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Abstract This is the first of a series of articles showing how 4 dimensionally co- 
variant analytical procedures developed in the context of General Relativity can be 
usefully adapted for application in a purely Newtonian framework where they provide 
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physical insights (e.g. concerning helicity currents) that are not so easy to obtain by 
the traditional approach based on a 3+1 space time decomposition. After an introduc- 
tory presentation of the relevant Milne spacetime structure and the associated Cartan 
connection, the essential principles are illustrated by application to the variational for- 
mulation of simple barotropic perfect fluid models. This variational treatment is then 
extended to conservative multiconstituent self gravitating fluid models of the more gen- 
eral kind that is needed for treating the effects of superfluidity in neutron stars. 



1. Introduction 

As a generalisation of previous work [1] on the special case of Landau's two- 
constituent superfluid model, using a 4-dimensionally covariant treatment of the kind 
pioneered by Peradzynski [2], this article presents a coherent fully covariant approach 
to the construction and application of Newtonian fluid models of the more general kind 
required in the context of neutron star phenomena in cases for which it is necessary to 
allow for independent motion of neutronic and protonic constituents. 

Whereas a simple perfect fluid model is sufficient for deriving the most basic features 
of neutron stars (such as the radius for a given mass and the oblateness for a given 
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angular momentum) models involving at least two independent constituents (of which 
one at least is superfluid) are needed to account for the details revealed by pulsar 
frequency observations. If quantitative accuracy is needed, the high mean densities of 
neutron stars require the use of a general relativistic treatment. In accordance with 
this requirement, in applications for which a simple perfect fluid model is sufficient, use 
of fully relativistic models has been standard practice since the outset of neutron star 
theory. However, when more elaborate models have been needed, most work has relied 
on less accurate Newtonian models, either because the relevant relativistic models had 
not been developed or because, even if available in principle, the relevant relativistic 
models were too difficult to apply in practice. 

When both non relativistic and relativistic versions are available, as is the case [1, 3] 
for multiconstituent superfluid models, the question of which is most appropriate for a 
given purpose depends not just on considerations of intrinsic accuracy or computational 
economy but also on questions of extrinsic compatibility with the relevant background 
framework. Thus for treating perturbations of a zeroth order global configuration de- 
scribed by a fully relativistic perfect fluid model, what will usually be most convenient 
is the employment [4, 5] of a two constituent fluid model that is also fully relativistic. 
However to deal with interactions with a solid crust described by a Newtonian elas- 
ticity model (since although appropriate relativistic elasticity models are available in 
principle [3] , their technical complexity has so far prevented them from being effectively 
applied in practice) it may be more practical [6, 7] to use a two constituent fluid model 
that is also non-relativistic. 

The purpose of this article is to show how to set up and apply a fully covariant 
formulation of the kinds of non-relativistic multiconstituent fluid dynamical models 
that are needed for such cases. Using a more traditional kind of formulation based on 
a preferred space reference frame (which complicates the treatment of effects such as 
helicity conservation, but facilitates the generalisation to allow for electromagnetism) a 
complementary development of the same class of Newtonian models has recently been 
provided by Prix [8]. 

The previous analysis on which the present work is based [1] was restricted to the 
Landau model which involves just a single massive particle constituent together with 
a second constituent, representing entropy, that is massless in the Newtonian limit. 
The more general analysis presented here covers cases (including the historic prototype 
of the Andreev Bashkin model [9] for a superfluid helium mixture) involving at least 
two independent constituents representing particles of which both kinds are massive, 
as in the particularly relevant example of the application to independently moving 
superfluid neutrons and protons. One of the side issues that will be dealt with here 
is the relationship between the effect of entrainment (whereby the relevant momenta 
deviate from the corresponding particle velocities) and the "effective masses" that have 
been defined in different ways in the published literature. 

As has already been pointed out by Peradzynski [2] , a noteworthy advantage (even 
for non-relativistic description) of using a 4-dimensionally covariant treatment, like that 
of the "canonical" approach developed here, is the possibility of exploiting Cartan type 
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methods of mathematical analysis, involving the use of antisymmetric differential forms 
for the construction, not just of vorticities, but also (as in the relativistic case [10]) of a 
more elaborate category of helicity currents. The technical intricacy of the construction 
procedure for these helicity currents is such that they would be very awkward to deal 
with (and are therefore usually ignored) in the conventional kind of approach based on 
a non-covariant frame formulation of Newtonian theory. 

When setting up a covariant description of Newtonian theory it is worthwhile to 
recall how, having already jettisoned the distinction between time and space in his spe- 
cial relativity theory, Einstein went on to jettison the distinction between spacetime 
and gravitation in his general relativity theory: in the latter theory the specification 
of the spacetime metric g^ automatically includes the specification of gravity via the 
corresponding Riemannian covariant differentiation operator V M which is determined 
by a connection with components f that are given by the well known Christoffel for- 
mula F u p = g U(T (g a (p, P ) — \gnp,a)-> using a comma for partial differentiation with respect 
to the (arbitrarily chosen) space time coordinates x M , and round brackets for index 
symmetrisation. However a distinction between spacetime and gravity can be made in 
the Newtonian analogue of Einstein's Riemannian structure, namely something that is 
known [11, 12] as a Newton Cartan structure - whose not so widely familiar principles 
will be recapitulated below - involving a non-Riemannian covariant differentiation op- 
erator D p that is determined by a Cartan connection with components uj " that are not 
obtainable from a Christoffel formula because the Newtonian structure does not spec- 
ify the non-degenerate spacetime metric that would be needed. In contrast with the 
inextricable case of general relativity, it is possible in the Newtonian case to extricate 
the specification of the gravitational field, as given by the (non-tensorial) Cartan con- 
nection components w", from that of the underlying spacetime manifold. Prior to any 
knowledge of the gravitational field (i.e. the Cartan connection) the only endowment 
of the underlying Newtonian spacetime manifold consists just of what is describable 
as a Milne structure, a concept that will be briefly recapitulated in the immediately 
following section. 

Like the Minkowski structure of special relativity theory, the Milne structure of New- 
tonian spacetime does not involve any free parameters or fields, i.e. it is intrinsically 
unique (modulo diffeomorphisms) . Nevertheless, despite its intrinsic simplicity the na- 
ture of this Newtonian spacetime structure is just sufficiently subtle to have prevented 
it from being properly understood until Milne's introduction of Newtonian cosmology 
theory a couple of decades after Einstein's introduction of general relativity -though not 
so long after Friedmann's foundation of the corresponding general relativistic cosmology 
theory - at about the same time as Cartan's epoch making work on the development 
of the appropriate mathematical machinery. Milne's breakthrough was based on the 
extrapolation to a global level of Einstein's earlier observation - originally in a Newto- 
nian context, as a guiding principle for the construction of the corresponding relativistic 
theory - of the equivalence at a local level between gravitation and acceleration. 
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2. Covariant description of Newtonian spacetime 

Whereas a fully covariant treatment is generally recognised to be indispensable for 
formulating general relativistic models, on the other hand, for their Newtonian ana- 
logues, the usual practice is to rely entirely on an 'Aristotelian' decomposition whereby 
space time is considered as a direct product of a flat Euclidean 3-space with a one 
dimensional Euclidean time line. Any such Aristotelian structure will be characterised 
by a corresponding class of Aristotelian coordinate systems, which consist of a set of 
ordinary Cartesian (orthonormal) space coordinates X 1 (i = 1,2,3) together with a 
Newtonian time coordinate t which is physically well defined modulo a choice of time 
origin. These coordinate systems are mapped onto each other by the transformations 
of a 7 parameter Aristotelian symmetry group, consisting of the product of the 6 pa- 
rameter group of Euclidean translations and rotations with the 1 parameter group of 
time translations. 

Whereas the time coordinate t of such a system is physically well defined (modulo an 
arbitrary adjustment of the origin) it has been generally recognised since the foundation 
of Newtonian theory that in a generic application there will be no uniquely preferred 
Aristotelian structure, but that the theory will be invariant with respect to a group of 
gauge transformations relating different Aristotelian product structures that all share 
the same constant time sections but that do not have the same sections of constant space 
position (as measured by fixed values of the Aristotelian space coordinates X 1 ). Any 
such gauge transformation will be specifiable by a mapping of one of the (7 parameter 
family of) sets of Aristotelian coordinates of a particular Aristotelian structure to a set 
belonging to another such structure. Since the flat constant time sections are preserved, 
any such transformation must be expressible just as a mapping X % i— > X\ t i— > t, for 
which the time transformation is trivial, t — t, and for which the new space coordinates 
X 1 are given by a time dependent Euclidean transformation. However not all kinds of 
time dependent Euclidean transformation are admissible. In particular time dependent 
rotations (belonging to what is known as the Coriolis group) are excluded from the 
status of gauge transformations because they change the physical comportment of the 
system (giving rise to what is known as the Coriolis effect). 

It turns out that the only admissible gauge transformations between different Aris- 
totelian structures are time dependent space translations as given by a transformation 
of the form X 1 i— > X 1 with 

X* = X i - c l , (1) 

for quantities c l that are constants in the sense that they have to be independent of 
the space coordinates X 1 , but that are arbitrarily variable as functions of time. For 
a long time it was generally believed that not all such time dependent translations 
were admissible, but only the three parameter set of Galilean transformations, meaning 
those that are linear so that the quantities c % can be taken to be given by expressions 
of the form c l = Ut in which the quantities b l are constants in the strong sense of being 
independent not just of the X 1 but also of t. A set of Aristotelian structures related by 
such linear transformations constitutes what may be described as a Galilei structure. 
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The important point that escaped everybody's notice until, after having been im- 
plicitly recognised in Einstein's "equivalence principle" , it was finally exploited for the 
foundation of Newtonian cosmology by Milne [14] is that Newtonian mechanics con- 
tains nothing to distinguish any particular preferred Galilei structure. The set of gauge 
transformations that are admissible in the sense that (unlike generic Coriolis trans- 
formations) they preserve the form of the physical laws of motion is not restricted to 
linear Galilean transformations, but includes generic transformations of the form (1) in 
which, while independent of the X\ the quantities c % are allowed to have an arbitrarily 
non-linear dependence on the time t. The complete set of all the Aristotelian product 
decompositions that can be obtained by such Milne gauge transformations constitutes 
what may be described as the Milne structure of space time. Thus the (physically 
unique) Milne structure consists of a family of (physically equivalent) Galilean struc- 
tures that are related to each other by accelerated space translation transformations of 
the form (1), while each member of the (infinite parameter set) of Galilei structures 
consists of a (3 parameter) set of Aristotelian structures that are related to each other 
by linear transformations of the form (1). 

Although the traditional employment of a particular choice of Aristotelian structure, 
and in particular of some corresponding set of Aristotelian coordinates {t,X 1 } is useful 
for many purposes, such as the exploitation of the flatness and distant parallelism of 
the preferred (constant time) 3-space sections, the advantages of this are not cost free. 
The price to be paid includes not just the well known obligation to verify Galilean 
(and Milne) invariance with respect to changes of the Aristotelian frame. Another, 
less well known, cost is the loss of access to the elegant and powerful mathematical 
methods based on tensors, and particularly on Cartan type differential forms, that 
become available when a fully covariant four dimensional framework is used. 

A convenient feature of the 3 dimensional constant time sections in an Aristotelian 
decomposition is the existence of a physically well defined - and flat - metric having 
components 7^ that are given simply by the unit matrix with respect to orthonormal 
Aristotelian space coordinates X 1 = 1,2,3). This metric can be used together 
with its contravariant inverse 7* J for raising and lowering of space indices, including 
those of the associated antisymmetric volume measure tensor e^k (whose non vanishing 

components are given by ±y / det{7} with the sign depending on whether the {i,j,k} 
is an even or odd permutation of {1,2, 3}) whose contravariant version will be charac- 
terised by the normalisation condition e^e^ = 3!. A similarly convenient feature of a 
general relativistic formulation, is the existence of a physically well defined - but not in 
general flat - space time metric giving components g^ v say with respect to coordinates 
(/i, v = 0,1,2,3), which determines a corresponding antisymmetric spacetime volume 
measure tensor with non vanishing components equal to —det{g}, and that can 
be used together with its contravariant inverse g^ v for raising and lowering of space 
time indices. Quite apart from the seductive flatness and parallelism properties of the 
preferred space sections in an Aristotelian decomposition, one of the reasons why fully 
covariant formulations of Newtonian dynamics are not as widely used as they deserve 
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to be is their lack of an analogous means of raising and lowering of space time indices. 
One of the purposes of this work is to show that, despite this handicap, Newtonian 
mechanics can nevertheless be set up without too much difficulty in a fully covariant 
formulation that makes it easy to exploit the technical advantages of freedom to use 
arbitrarily curved space time coordinates x M . 

The first step is to obtain the fundamental spacetime tensor fields that are available, 
in lieu of the relativistic metric tensor, for the characterization of Newtonian space 
time. To start with, for any fixed value of the preferred Newtonian time coordinate, t, 
the embedding mapping X 1 \— > x M of the corresponding Aristotelian space section will 
determine physically well defined contravariant spacetime tensor fields with components 
given by 

T = l ij xWi = e^x^Xk , (2) 

using a comma to indicate partial differentiation (in this case with respect to the space 
coordinates X 1 ). However there is no unambiguously preferred (Galilei invariant) pre- 
scription for lowering the indices to obtain corresponding contravariant versions 7^ and 
^livp because the former is degenerate (7^ has matrix rank 3 not 4) so it does not have 
a well defined inverse. 

The foregoing consideration means that in Newtonian theory 4-dimensional tensor 
indices will in general have an irrevocably covariant or contravariant nature. The most 
basic example is that of the preferred covariant unit vector £ M that is obtained simply 
as the gradient of the preferred Newtonian time coordinate, i.e. = t tfM , and that is 
a null eigencovector of the degenerate preferred contravariant space metric, and also 
orthogonal to £^ p , i.e. 

7"% = 0, eT% = Q. (3) 

An exception, having both a covariant and a contravariant version, is that of the anti- 
symmetric spacetime volume measure tensor £ Mi , pcr , which has a physically well defined 
normalisation - despite the lack of a non degenerate space time metric in the Newtonian 
framework - that is specified by the relation 

and for which a corresponding contravariant version is unambiguously definable by a 
normalisation condition of the same form, 

£p,vpo-£^ ^ 4! , (5) 

as holds in general relativity, which means that it will satisfy 

eT^ta = e^ p . (6) 

The fields 7^ and £ M (the degenerate residue representing all of the algebraic struc- 
ture that remains from the relativistic spacetime metric in the Newtonian limit) con- 
stitute what may be termed a Coriolis structure, since it contains nothing that distin- 
guishes rotating from non rotating frames. To incorporate this distinction, and thereby 
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complete the covariant specification of Newtonian space time, we must consider what to 
use in place of the Riemannian connection and the associated covariant differentiation 
operator that in general relativity is uniquely specified by the spacetime metric g^ v . 
In the Newtonian case, a corresponding physically well defined but non-Riemannian 
connection u " and associated covariant differentiation operator will be provided 
by the Newton-Cartan structure that is determined by the gravitational field in the 
manner to be described at the end of the next section. However the only features of 
the Newtonian space time background that are well defined a priori (in the absence of 
information about the gravitational field that specifies the Newton -Cartan structure) 
are those provided by the Milne structure that was described at the beginning of this 
section. 

The Milne structure by itself (without reference to the gravitational field) does not 
provide enough information for the unambiguous specification of a connection. What 
it is capable of providing is a connection that is gauge dependent. Specifically for 
each choice of gauge, i.e. for each choice of a particular Aristotelian product structure, 
there will be a corresponding natural connection with components f defined by the 
condition that they simply vanish when evaluated with respect to a corresponding 
system of Aristotelian coordinates {t,X 1 } (though not of course if the orthonormal 
space coordinates X 1 were replaced by coordinates of some other, e.g. spherical, kind). 
This means that with respect to coordinates of this particular type (but not for a system 
of some other, e.g. spherical, kind) the corresponding covariant differentiation operator 
V M will be identifiable with the simple partial differentiation operator d^. 

The gauge dependent connection r " that is defined in this way has two convenient 
properties. Firstly, since it is identifiable with partial differentiation in the relevant 
Aristotelian coordinate system, it is clear that since d^d u = d u d^ the covariant derivative 
will automatically inherit the analogous commutation property 

V M V, = V,V M , (7) 

i.e. the connection F v has the property that (unlike the curved Newton - Cartan 
connection described below) it is flat. The other convenient property is that although 
(again unlike the gravitational field dependent Newton-Cartan connection oj"^ it is 
gauge dependent, its gauge dependence is of a rather weak kind. Since connection 
components are unaffected by linear transformations, it follows that f will not be 
affected by gauge transformations of the restricted Galilean type. Choosing a partic- 
ular Galilean structure (i.e. a linearly related subclass of Aristotelian structures) is 
equivalent to choosing a particular connection of this flat type. It will be shown below 
how such a connection is affected by Milne gauge transformations of the more general 
accelerated type that relate distinct Galilean structures. 

According to the preceeding definition, the Milne structure is an equivalence class 
of Aristotelian (direct product of time and flat space) structures that are related to 
each other by a gauge group consisting of time dependent space translations. In the 
modern (post Cartan) technical language of fibre bundle theory this definition can be 
succinctly reformulated as follows. 
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In formal mathematical terms, the Milne structure of Newtonian spacetime is for- 
mally describable as bundle of 3-dimensional Euclidean space fibres (each characterized 
by its own flat metric 7^) over a base consisting of a line endowed with a physically 
preferred time measure (specified by the coordinate t) for which the relevant (Abelian) 
bundle group consists just of the 3 - parameter set of Euclidean space translations (but 
not rotations), which are expressible in terms of Cartesian coordinate X % on the fibre 
simply by transformations of the form (1). (The exclusion of Euclidean rotations from 
the bundle group is what distinguishes the Milne structure from the more primitive 
Coriolis structure given just by the specification of the fields and 7^.) 

Any particular section of this Milne bundle, i.e. a representation as a direct product 
of the base times the fibre, will be interpretable as a particular choice of an Aristotelian 
structure. In such a direct product structure, the preferred time coordinate t on the 
base and a choice of Cartesian coordinates X 1 on a space section (i.e. coordinates such 
that metric components 7^ form a unit matrix) will determine a corresponding set of 
Aristotelian spacetime coordinates x M according to the obvious specification x° = t, 
x % = X 1 . A corresponding connection is thereby definable as the one with respect 
to which covariant differentiation reduces to partial differentiation, i.e. the one for 
which V M = dfj, and T p v = 0, in these particular coordinates. (This connection could 
be interpreted as the Riemannian connection provided by a flat Unruh type spacetime 
metric of the form ds 2 = 7^ dX l dX j — C 2 dt 2 in which C could be any arbitrarily chosen 
constant speed, which might be that of light, but which in applications to perturbations 
in an asymptotically uniform fluid background could more usefully [13] be taken to be 
the relevant sound speed.) 

3. Galilean and Milne type gauge transformations 

Although the dynamical equations of ordinary Maxwellian electromagnetic theory 
are expressible in terms just of gauge independent quantities such as the electric and 
magnetic fields, it is useful for many purposes, and indispensible for a variational formu- 
lation, to employ various gauge dependent entities (starting with the vector potential 
Ap) . In the present context (of ordinary Newtonian dynamical theory) it is analogously 
useful for many purposes, and indispensible for a variational formulation, to employ en- 
tities whose specification depends on a particular choice of gauge, where in this context 
a "choice of gauge" is to be interpreted as meaning a particular choice of Aristotelian 
structure within the large equivalence class of Aristotelian structures that collectively 
constitute the Milne structure described above. 

In terms of the Aristotelian coordinates {t, X 1 }, the transformation to the analo- 
gous coordinates for any alternative bundle section, i.e. any alternative Aristotelian 
structure, will be expressible by a transformation in which the base coordinate is held 
fixed, i.e. 1 1— > t while the Cartesian fibre coordinates X 1 are transformed according to 
a relation of the form (1) in which the translation vector d is given as an arbitrarily 
variable function of the base coordinate t. It is evident that the connection specified by 
the new gauge will be identical with the one specified by the old section, i.e. we shall 
have f ' " = V u p and hence V M = V M , so long as the transformation (1) is linear, i.e. so 
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long as the translation is of the Galilean form characterised by the condition that the 
components 

V = ^, 6o=0, (8) 

of the boost velocity vector of the transformation should be fixed, independently of t. 
However the connection will not be preserved if there is a non-vanishing value for the 
corresponding relative acceleration vector a M as given by 

a l = -^, a = . (9) 
dt 

It can be verified that for the new section obtained by a generic (i.e. accelerated) 
Milne gauge transformation the corresponding new connection will be related to the 
original one by a relation of the simple but non trivial form \— > F " p with 

V P = r/, - , (io) 

where a u is the relevant transformation generator, which can be any vector that is 
spacelike and spacially uniform (i.e. purely time dependent) in the sense that 

= , 7 ^VX = . (11) 

This transformation law has the noteworthy feature of preserving the trace of the con- 
nection, i.e. it gives 

f/„ = t; v . (12) 

It follows that if n M is a physically well defined current 4-vector of the kind to be 
discussed in the next section then the divergence given by the expression V u n u will also 
be physically well defined as a gauge independent scalar field (which will vanish in the 
particular case of a current that is conserved). 

It will be convenient for future reference to introduce a scalar boost potential function 
/3 that is defined, modulo an arbitrary time dependent constant of integration f3{0} by, 

If = rfV v p , ^VJf = , (13) 

so that in the original Aristotelian coordinate system it will be given explicitly by an 
expression of the form 

P = P{0} + ^VX\ 7 ^V„/3{0} = 0. (14) 

It then follows that the relative acceleration vector will be given by 

a» = e u V u lf = 7^V„ a, a = e u V u (3 , (15) 

where e M is the relevant "ether" 4-velocity vector, i.e. the unit time lapse vector of the 
Aristotelian rest frame, whose components with respect to the corresponding coordi- 
nates {t, X % } will be given by e° = 1, e l = 0, so that it will satisfy the conditions 

e% = l, VX = 0. (16) 



9 



Whichever bundle section may have been used to specify it in the first place, the 
connection will automatically be such as to preserve the fundamental tensors 7^ and 
of the (Coriolis) spacetime structure, i.e. it will satisfy 

V m7 " p = 0, V^ = 0, (17) 

and hence also 

V^ paT = . (18) 

Conversely the specification of any particular connection that satisfies these preserva- 
tion conditions will characterize what may be described as the corresponding Galilei 
structure, which can be conceived as an equivalence class of Aristotelian (i.e. direct 
product) structures related by linear gauge transformations of the form (1) with van- 
ishing value of the relative acceleration vector defined by (9). 

For centuries after Newton's original development of his theory it was taken for 
granted that the relevant equations of motion singled out a preferred Galilei structure 
with a corresponding Galilean transformation group with respect to which their form 
remained covariant. What Milne realised [14] was that except in the case of a localised 
system in an asymptotically empty background (such as the solar system example to 
which the early successes of Newton's theory were restricted) the equivalence princi- 
ple prevents the prescription of any natural rule for preferring some particular Galilean 
structure rather than another. Thus, as a consequence of the applicability of the equiva- 
lence principle, it transpires that the relevant equations of motion are covariant not just 
with respect to the Galilei group but also with respect to the larger Milne group, which 
relates distinct Galilei structures by transformations characterized by non-vanishing 
values of the relative acceleration vector a M . 



4. Gravity, particle dynamics, and the Newton - Cartan connection 

The way the foregoing considerations apply to the most basic of the Newtonian 
equations of motion, namely the equation of motion for a free particle in some given 
gravitational field, is as follows. Having specified the worldline of the particle giving 
the relevant spacetime coordinates as functions of the Newtonian time t, one can 
go on to define the corresponding 4-velocity vector defined as the time parametrised 
tangent vector given by 

« = (19) 

of which only three components are actually independent since the definition automat- 
ically ensures that it satisfies the unit normalisation condition 

u% = 1 . (20) 

The equation of motion will then be expressible covariantly as 

u v V v u» = g», (21) 



10 



where is the relevant gravitational field vector, which must be strictly spacelike to 
avoid inconsistency with (20), and which, more specifically, is required to be derivable 
as the space gradient of a Newtonian potential 0, i.e. 

W = > Sf = -7^,0 . (22) 

Although it was traditionally expressed in a non-covariant mathematical form, the 
relation (21) was recognised from the outset to be physically invariant with respect to 
Galilean transformations, i.e. the linear transformations that preserve the connection 
involved in the covariant differentiation operator V^. The crucial point that eluded 
Newton and everyone else before the time of Einstein, Friedmann, and Milne is that 
is akin to the scalar potential 4> (and to Maxwell's covector potential as contrasted 
with the tensorially well defined electromagnetic field F^ v ) in that it cannot be con- 
sidered to be an absolutely well defined locally measurable vector field but should be 
recognised to be gauge dependent. It can be seen from (10) that the relation (21) is 
in fact invariant not just with respect to Galilean transformations but to generic Milne 
transformations as characterised by a non vanishing acceleration vector in (9), pro- 
vided it is understood that g^ undergoes a corresponding Milne gauge transformation 
of the simple form g M i— > g^ with 

= 0" - a" . (23) 

This evidently entails the requirement that the Newtonian potential should transform 
according to a law of the form <fi i— > </> with 

4> = <t> + a , (24) 

in which scalar field a will be given by (15) as the ether (i.e. Aristotelian) frame 
time derivative of the boost potential j3 . The freedom to adjust the specification (14) 
of the latter by freely choosing the time dependence of the value /?{0} of the boost 
potential at the Aristotelian coordinate origin corresponds to the calibration freedom 
in the specification of by (22). Thus, even in a fixed Aristotelian frame, the potential 
<fi will be subject to trivial gauge transformations consisting of the addition of a purely 
time dependent quantity that is identifiable simply as the time derivative of /?{0}. 

The idea of the Newton Cartan formulation [11, 12] is to replace the gauge depen- 
dent differential operator V M by a corresponding gauge covariant differential operator 

that is specified by replacing the flat connection by a gravitationally modified 
connection a; " say that is given by 
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Y; p - t,g% . (25) 



Using (23) in conjunction with (10), it can be verified that this modified connection 

3 erty w/ p i-> u» p with 

<Vp = < p • (26) 



has the desired gauge invariance property u^" i— > u " with 
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This makes it possible to rewrite the dynamical equation (21) in the manifestly gauge 
invariant (as well as coordinate covariant) form 

u u D v u» = . (27) 

While facilitating the comparison with general relativity, this gauge covariant differen- 
tiation operator has the disadvantage of lacking the convenient flatness property 
(7) of the gauge dependent alternative V M . For example if we consider not just a single 
particle trajectory but a fluid flow characterized by a velocity 4- vector that is defined 
as a field over spacetime then we shall have 

W=|V/, (28) 

using square brackets to indicate index antisymmetrisation, where R^ a is the Newton- 
Cartan curvature tensor which can easily be seen to be given by the expression 

R fl f a = 2t a t [ „V l/] g». (29) 

Although it is not immediately manifest from this formula, it can be easily verified 
using (11) and (23) that this curvature tensor is indeed invariant under the gauge 
transformation (10), i.e. it satisfies R fll f a i— > R„ v p a with 

K p o = v* • ( 3 °) 

It can be seen that the corresponding (again Milne gauge independent) Ricci type 
curvature trace tensor is proportional to the Laplacian of the gravitational potential, 
having only a single independent component that is specified by the formula 

R^ = K*\ = V-7 PCT V P V CT 0. (31) 



5. Action and the 4-momentum covector 

Although, as has just been shown, the dynamical equation (27) is gauge independent, 
its derivation from an action principle requires the use of a gauge dependent momentum 
covector 7^ that is in many ways analogous to the electromagnetic gauge potential 
that is needed for the variational formulation of Maxwell's equation. The original 
variational formulation of the Newtonian dynamical equation by Laplace was given by 
an action integral 

1 = J Ldt, (32) 

in which the Lagrangian scalar is taken to be the difference between the kinetic and 
potential energies in the well known form 

L = \mv 2 — mcp , (33) 



12 



where m is a constant mass parameter, and v is the magnitude of the 3 velocity vector 
as specified with respect to some chosen Aristotelian reference system, in terms of 
the corresponding ether vector e^, by 

^ = u » - e » , v% = , (34) 

so that its Aristotelian components will be given by 

A Vi 

v" = , v l = . (35) 
at 

and its magnitude by 

v 2 = 7ijU V . (36) 
This action can be rewritten in the more elegantly covariant form 

1 = J tt^ dx^ 1 , (37) 

which is evidently equivalent to taking 

L = ttX , (38) 

by defining the appropriate gauge dependent 4-momentum covector as follows. In terms 
of the Aristotelian coordinate system x° — t, x l — X 1 characterizing the chosen gauge, in 
which we shall evidently have 

u° = 1 , = (39) 

the appropriate 4-momentum covector will be given by the prescription 

7r = -£, TTi = m^ijV 1 , (40) 
in which the total particle energy is given by the usual formula 

£ = |mv 2 + mcf) . (41) 

Instead of explicit reliance on the Aristotelian coordinate system, we can use the cor- 
responding ether vector e M as introduced in (16), i.e. the unit 4- velocity vector of the 
Aristotelian rest frame, with respect to which its components will be given simply by 

e° = 1 , e = , (42) 

for the purpose of obtaining a covariant expression for 7r M as follows. To start with, 
we use the Kronecker unit tensor 5^ u to construct the corresponding Aristotelian space 
projection tensor according to the specification 

7<* = n " e% • (43) 
13 



We then use the defining relations 

l„e v = , 7^7^ = , (44) 

to characterize the covariant tensor 7^ obtained via the Aristotelian coordinate map- 
ping from the space metric 7^ in the constant time 3-sections. We can then define the 
kinetic 4-momentum vector by 

Pn = mv^ - \mv 2 t il , (45) 

where 

*V = lw uV i v 2 = v ^ = l^u v (46) 

so that we have 

1 .2 
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mv = p u u u = —p v e v . (47) 



Like the the kinetic momentum has only three independent components, being sub- 
ject to a constraint that in this case (unlike that of -u M ) is ether frame dependent, having 
the form 

Tv^Vv = -2me% . (48) 
In terms of this kinetic contribution we finally obtain the expression 

^=Pn-m4t IM , (49) 

for the complete momentum covector. 

6. Finite and infinitesimal gauge transformation rules 

Since the boost transformation law for the Aristotelian ether vector evidently takes 
the form e M 1— > e M with 

e» = + If , (50) 
it can be seen that while, the degenerate contravariant metric tensor is gauge invariant, 

*T = 7"" , (51) 

the corresponding mixed projection tensor (43) will undergo a transformation given by 

it = -fu~ tuP , (52) 

and the corresponding degenerate covariant metric tensor will be governed by the less 
trivial transformation rule 

l»v = l»v - 1t{plv) P b p + , (53) 

in which the boost magnitude b is naturally defined by 

ft 2 = 7^6". (54) 
14 



Thus while the contravariant form of the relative velocity (34) obeys the simple Galilean 
transformation rule 

v» = - If , (55) 

the corresponding covariant vector 

= l^v v (56) 
transforms according to the less simple rule 

v lt = v ll - j^b u + t^b 2 - b v v v ) , (57) 
while for the squared velocity 

v 2 = u% = -y^vPu" , (58) 

one obtains 

v 2 = v 2 - 2b^ + b 2 . (59) 

When applied to the kinetic momentum covector (45) the foregoing formulae provide 
the gauge transformation rule h- > with 

Pn = Pp- fniiivb v + |m&% , (60) 

and thus the transformation rule for the complete momentum covector (49) can be seen 
to be expressible in terms of the boost potential (3 by 

= TT/i rnV^p + \m b\ , (61) 

so that, in particular, the corresponding transformation law for its energy component 
(41) will be expressible as 

£ = S - Pv b v + m(\b 2 + a) , (62) 

while finally for the Lagrangian scalar (38) we obtain L i— > L with 

L — L — mu"V M /3 + \m b 2 . (63) 

It is apparent from the form of these transformation rules that it will be convenient 
to work with a recalibrated boost potential, f3 — f3 — \ J b 2 dt, that will be characterized 
by 

Vj = V u f3 - \b 2 t u . (64) 

Since the difference between j3 and is a function only of the cosmological time, we 
can just as well use the latter as the former in the characterization (13), which can be 
rewritten as 

If = rfV v $ , 7"" W = , (65) 
15 



but in terms of the recalibrated boost potential the expression (15) for the acceleration 
will acquire the slightly modified form 

a » = -fV v a , a = efV^ + \b 2 . (66) 

The corresponding modification of the formula (62) for the transformation of the energy 
will be given by 

£ = £- Pv b v + m(e u + b")Vj . (67) 

However it is for the 4-momentum covector and the Lagrangian that the advantage 
of the modified boost potential becomes apparent, since it allows (61) and (62) to be 
rewritten more simply as 

= tt m - rriW^(3 , (68) 

and 

L — L — mu^V^p . (69) 

The effect of the gauge transformation on the action integral (32) can thus be seen 
to be given by X i— > X with 

X = l-m[(3], (70) 

using large square brackets to denote the total change in the boost potential /3 along 
the worldline segment under consideration. Since the gauge adjustment term in (70) 
will not be affected by any purely local variation of the worldline (local meaning that 
it is non vanishing only on a confined subsegment lying entirely inside the extended 
worldline segment under consideration) it is obvious that it will have no effect as far 
as the application of the variational principle is concerned. The observation that the 
gauge transformation changes the action only by an amount that is constant in the 
sense of being independent of local worldline variations evidently accounts for the in- 
variance with respect to these (linear Galilean or accelerated Milne) transformations of 
the ensuing system of dynamical equations. 

The preceeding example is one of many in which it is simpler not to work with finite 
(Galilean or Milne) transformations as we have been doing so far, but with linearized 
infinitesimal transformations. For a given gauge transformation generated by a given 
boost potential /3, whereby a generic quantity, q say, is subject to a mapping q i— > q, 
the corresponding infinitesimal gauge transformation g i— > dg is defined by a routine 
two step procedure as follows. The first step is to construct a homotopic interpolation 
by a one parameter family of gauge transformations q h- > q{e} with q{0} = q and 
<?{!} — 0.1 f° r which q{e} is given, for intermediate values of the homotopy parameter 
e, by an interpolating boost potential (3{e} = e/3. The corresponding infinitesimal 
transformation is then obtained by taking the limit, as e — > 0, of the derivative with 
respect to e, so that we have 

a « = !"S5l(<*»- < 71 > 

This differential operation will merely restore the value we started with (as recovered 
by setting e = 0) for quantities whose dependence on the transformation amplitude is 
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homogeneously linear as is the case for the diverse derivatives of the boost function, 
which will be characterized by 



da = a, d> = ^, da" = a", (72) 

and also of course for the original boost function (3 itself, though not for the modified 
boost function (3 for which we obtain a convenient simplification, 

dp = dfi = P, (73) 

which means that at the differential level the modification (3 is redundant. In a similar 
manner, going to the differential level provides no great simplification for quantities 
whose gauge transformation depends just linearly on the boost amplitude, as is the 
case for ether vector e M , and the 3- velocity vector v^, for which we obtain 

de" = If , d> = -If , (74) 

as well as for the less trivial cases of the mixed projection tensor 7^, and the connection 
ly^, for which, from (52) and (10), we obtain 

d 7 £ = -t v \f , dry; = -t, a % , (75) 

but for quantities with more complicated gauge dependence the differential level is much 
more convenient. Noteworthy examples are the covariant space metric tensor 7^ and 
the kinetic momentum vector for which the formulae (53) and (60) reduce simply to 

cb> = -^t{ l xlv) P b p , dp M = -my^V , (76) 

and particularly the Lagrangian L, and the complete momentum covector 7r M , for which 
we simply obtain 

dL = -mu u V u (3 , d7r M = -mV^ . (77) 

(The feature of transforming just by the addition of the gradient of a scalar is a property 
that the 4-momentum covector 717 of Newtonian dynamics shares with the covector 
potential of Maxwellian electromagnetism, a relationship that is essential for the 
amalgamation of gravity and electromagnetism in the corresponding general relativistic 
theory [3].) 

The corresponding differential version of the formula (70) for the action integral is 
given by the expression 



dX = — m 



(3 , (78) 



whose evident path independence makes it obvious that the ensuing dynamical theory 
will have to be gauge invariant, as shown above by the existence of the manifestly 
covariant formulation (27) of the equations of motion. 
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7. Covariant fluid current variation formulae. 

The simple Lagrangian worldline variation principle for a single particle that was 
discussed in the preceeding sections can be generalised in an obviously natural way to 
a flow line variation principle for a fluid system. 

Before proceeding, it is to be recalled that in the absence of any prescribed spacetime 
structure the usual description of currents in terms of vector fields will not be available, 
but it will still be possible to use the more fundamental Cartan type description whereby 
a current is represented as a 3-form, i.e. an antisymmetric covariant tensor with 4- 
dimensional components Np Up whose surface integral 

N=^J Np Up d 3 x^ (79) 

determines the total number flux over a 3-surface with tangent element d 3 x pup = 
3\d {1) x^ d {2) x u d {:i) x p } generated by infinitesimal displacements d (i) x p ( i— 1,2,3). In such 
a description, the condition for the conservation of the number flux is the vanishing of 
its exterior derivative as defined by (d A N)p Upa = 4d[pN upa ]. A convenient feature of 
such exterior differentiation is that it makes no difference if the partial differentiation 
operator dp is replaced by a tensorially covariant differentiation operator V M (or by Car- 
tan's gauge covariant differentiation operator Dp) since due to the antisymmetrisation 
all the (symmetric) connection components will cancel out. 

The more fundamental three index description of the current will of course be repla- 
cable by a more compact description involving a single contravariant index whenever 
the spacetime background is endowed with a canonical antisymmetric measure tensor 
Ep Vpa and a corresponding contravariant alternating tensor e pvpu ^ which, as seen above, 
will be the case both in relativistic and Newtonian spacetime. The required current 
vector, with Aristotelian rest frame components 

n° = n , n l = nv i , (80) 

where n is the ordinary (scalar) particle number density, will then be given by the 
duality relation 

n , = l £ ^N vpa , (81) 

which can be inverted to provide the expression 

Np Vp = Ep Upa rf . (82) 

For any covariant symmetric connection compatible with the measure preservation con- 
dition (18), the corresponding covariant differentiation operator will determine a diver- 
gence that will just be the dual of the exterior derivative operator. Thus independently 
of the choice of gauge one obtains the equivalent expressions 

VX = DX = T,^ upa (d A N)p Upa = ^-e^N^p , (83) 
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for the particle rate, which will vanish for a current that is conserved. 

In order to apply the variational principle, we need to evaluate the variation of the 
current that will result from transport by the action of an infinitesimal displacement 
x ii ^ x n _|_ As a general principle [3] the fixed point (Eulerian) variation of any field 
will be given generally by its comoving (Lagrangian) variation minus the relevant Lie 
derivative. Since the specification of the covariant representation of the current does 
not depend on any background structure its comoving variation will simply vanish, so 
its fixed point variation will just be given by the formula 

SN^ = -££AW> , (84) 

in which the Lie derivative is given by the standard formula 

&N^ p = i a W a N^ p + 3N a[llu V p] C • (85) 

As for any Lie derivative formula (and as in exterior differentiation) it makes no 
difference if the partial differentiation operator <9 M is replaced by a tensorially covariant 
differentiation operator (or by Cartan's gauge covariant differentiation D M ) since due 
to the antisymmetrisation all the (symmetric) connection components will cancel out. 
Thus by taking the dual of the formula (85) we obtain the useful theorem that for any 
covariant symmetric connection compatible with the measure preservation condition 
(18), the corresponding covariant differentiation operator can be used to express the 
fixed point (Eulerian) current variation produced by the displacement vector field ^ in 
the form 

Sn" = n u V v ^ - e^.n" - ra"V„f " . (86) 

This result establishes the validity in a Newtonian framework of a formula that has 
long been in regular use in a relativistic context, where it was originally derived by a 
rather different line of reasonning [3] that depended on the Riemannian specification of 
the covariant differentiation in terms of the non-degenerate space time metric that is 
no longer available in the Newtonian case. It is clear from this present approach that 
the formula (86) will remain valid if Cartan's gauge covariant derivative operator D p is 
substituted in place of the flat but gauge dependent derivative operator V M . 

A useful corollary of (86) is the corresponding formula for the variation of the current 
divergence, which takes the simple form 

8& v n v ) = -V^V v n v ). (87) 

It is immediately apparent from this that if the original current is conserved, i.e. if 
V u n u vanishes, then the displaced current will have the same conservation property. 

The formal identity of the relativistic and Newtonian variation formulae will be lost 
if we make a decomposition of the usual form 

n 11 = nu /1 , n = n% , (88) 

in which n is the ordinary particle number density scalar and is the 4-velocity of 
the flow as characterized by the unit normalisation condition (20). By contracting (86) 
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with it can seen that the variation law for the particle number density will be given 
by 

5n = tXV„e ~ V„(rO , (89) 

which has a different form from its relativistic analogue [3] due to its dependence on 
the preferred time basis vector (instead of the non degenerate spacetime metric 
that plays the corresponding role in the relativistic version). The same remark applies 
to the corresponding variation of the 4- velocity of the flow, which can be seen from (86) 
and (89) to be given by 

Su" = u u V„e - fvy - v! i u v t p Vu£ p ■ (90) 



8. Action principle for simple perfect fluid. 

The natural way to extend the single particle action principle discussed above to a 
corresponding fluid action principle is to base the latter on the space time integral of a 
Lagrangian density A that will be given by a decomposition of the form 

A = A POt + A kin + A int , (91) 

in which the first two (gauge dependent) terms are given by the product of the relevant 
particle number density n with the corresponding single particle contributions, while 
the extra term A. nt is given simply by 

A t = ~U int , (92) 

where U. nt is the internal compression energy, which will be a (naturally gauge indepen- 
dent) function just of the particle number density n. Specifically, the external potential 
energy contribution will be given by usual Newtonian formula 

A pot = -nmcf) , (93) 

which may be rewritten in covariant form as 

A pot = -0p%, (94) 

where the mass density current is defined by 

p M = pu^ = mn^ , p = nm , (95) 

while in accordance with (47) the kinetic contribution will be given by 

A kin = n vX = n% , (96) 

in terms of the (gauge dependent) kinetic 4- momentum covector defined by (45). The 
extra (gauge independent) contribution (92) representing the negative of the internal 
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compression energy density will determine a corresponding (gauge independent) chem- 
ical potential \ by a variation formula of the standard form 

5U int = X Sn . (97) 

In terms of this chemical potential function, the relevant perfect fluid pressure function 
P (which is also gauge independent) will be given by the well known formula 

P = n X - U int . (98) 

Since, as a consequence of the restriction (48), the variation of the kinetic momentum 
is automatically constrained to satisfy the identity 

u^p^ = , (99) 

it follows that the generic variation of the purely kinetic contribution to the Lagrangian 
will be given simply by 

5A kin =p,5n», (100) 

The variation of the combination (91) will therefore be given by an expression of the 
canonical form 

8 A — r K jl bn lx — p5(f) , (101) 
in which the total 4-momentum is given by an expression of the form 

7i> = - m</> £ M , (102) 

which differs from the corresponding free particle momentum formula (49) by the re- 
placement of the purely kinetic contribution p^ by a total material 4-momentum cov- 
ector \Xp, that is defined by 

Vn=Pn- Xt^ ■ (103) 
This material momentum covector is alternatively definable by the variation formula 

M mat = /v^, (104) 

where the total material (non gravitational) contribution to the Lagrangian density is 
defined by 

A = A. . + A. . . (105) 

mat kin int V / 

It can now be seen that the complete Lagrangian (91) will be elegantly expressible 
in terms of the 4-momentum covector ir^ and the pressure function P as 

A = n% + P. (106) 

It is to be remarked that while the pressure term in (106) is gauge invariant, the first 
term is not. However as the formula (102) for the fluid particle 4-momentum covector 
7r^ differs from its single particle analogue (49) only by the gauge independent term 
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proportional to \ in (103), it can be seen its variation dV^ under the action of an 
infinitesimal gauge transformation will be given by the same simple formula (77) as in 
the single particle case. It follows that the corresponding infinitesimal gauge variation 
of the Lagrangian density (106) will be given simply by 

dA = -p u V u p, (107) 

If it is taken for granted that the fluid obeys the ordinary Newtonian mass conservation 
law 

Vuff = 0, (108) 
the infinitesimal gauge variation will be expressible as a pure divergence in the form 

dA = -V„(/?,0 . (109) 

This means that the gauge change will have no effect on a localised variation of the 
space-time volume integral of the action density, so the dynamical equations given by 
the action principle will automatically be gauge independent. 

For the actual evaluation of the variation of the action, the work of the preceeding 
section provides all the elements that are needed. It can be seen from (101) and (86) 
that when the fluid flow is subjected to the action of an infinitesimal displacement 
vector field the resulting variation of the Lagrangian density will be given by the 
formula 

5A = v^^r 1 ) - - ps<f> , (no) 

in which the covector f ^ is interpretable as the 4-force density acting on the fluid, 
excluding the gravitational contribution which is already taken into account within the 
formalism. This force density can be seen to be given by the prescription 

f li = 2n v V [v Tr lt] +Tr lt V v n v , (111) 

i.e. it is constructed from the current vector and the corresponding momentum covector 
by contraction with (exterior) derivative plus derivative of contraction. 

The postulate that the mass parameter m should be constant means that the mass 
conservation law (108) will be equivalent to the particle conservation law given, in 
accordance with (83) by 

V„ra" = 0, (112) 

a result that is alternatively derivable, wherever the product n^n^ is non zero, as a 
consequence of the variational requirement that the force density (111) should vanish. 
Subject to (112), the expression for the force density will evidently reduce to the simple 
form 

fp = n"w vll , (113) 

in which the relevant generalized 4-vorticity 2-form is defined, as the exterior derivative 
of the 4-momentum covector, by 

07^ = 2^7^]. (114) 
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This vorticity 2-form w pv is generalized in the sense that it automatically includes 
allowance for gravity whose effect can be separated out in the decomposition 

tZfw = V + 2mt [fj y u ](fi . (115) 

in which is the ordinary material vorticity tensor defined by 

w pu = 2V[^] . (116) 

The adjustment allowing for gravitation affects only the time components, so both the 
complete and the material vorticity give the same purely spacelike 3-velocity vector, 
which is expressible independently of (p as 

w » = \ £ ^ Wvp = \s^w vp , (117) 

and which is related to the purely kinematic local angular velocity vector by the 
proportionality relation 

w» = 2mw^ , up = \e^ p V u v p . (118) 

A related quantity that is easy to analyse in the covariant formalism we are using 
here, but much more awkward to treat using the traditional 3+1 spacetime decompo- 
sition - so much so that its role in Newtonian fluid dynamics, was not recognised until 
the relatively recent work of Moreau [15] and Moffat [16] (a century after the pionnering 
analysis of vorticity by workers of Kelvin's generation) is that of helicity. In that work 
(and in its more recent non-barotropic generalisation [17]) the helicity was introduced 
as a scalar density that was constructed as the three dimension scalar product of the 
velocity vector v % and the vorticity vector w l . On the basis of experience [3] with the 
relativistic case, it is evident that in the 4-dimensionally covariant formalism we are 
using here, the helicity will most naturally be definable [2, 1] as a vectorial current 
rf that is proportional to the dual of the exterior product of the energy momentum 
covector ix p with the corresponding generalized vorticity two form, w pvi namely 

rf = E^7r u V p n a = \e^°* v w pc . (119) 

The time component of this quantity can be seen from (6) to be proportional to the 
Moreau-Moffat helicity scalar, with a negative coefficient (in the sign convention we are 
using, for which the sign of the measure component e 0123 is taken to be positive, so that 
of e 0123 will be negative) that is given by —2m 2 , i.e. 

rf tp = -w^ty p = -2m 2 Jvi . (120) 

It immediately follows from the Eulerian dynamical equation, whose variational deriva- 
tion will be described below, that this helicity current rf will be conserved in the simple 
sense that its 4-divergence, 

V/ = V^ pffl (121) 
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will turn out to vanish, a property that is both obvious and easy to express in the 4- 
dimensional approach used here, but not so trivial, either to derive or even to present, 
within the (Latin as opposed to Greek index) framework of the traditional 3+1 formal- 
ism (see appendix). 

The purport of the variational principle is that the spacetime volume integral of 
5A should vanish for any displacement £ M with bounded support (i.e. that vanishes 
outside some bounded spacetime region). Since, by Green's theorem, the divergence 
contribution in (110) will make no contribution to the variational integral, the principle 
reduces to the requirement that the 4-force density should vanish, 

U = 0. (122) 

Subject to (112) this equation will reduce to the form 

n v w vli = 0, (123) 

in which the complete vorticity 2-form can be expressed as 

t*V = 2V I^] + 2t l^u] (X + . (124) 

An immediate consequence of (123) is that (since its components form an antisymmetric 
matrix with a zero eigenvalue eigenvector, namely n^) the vorticity 2-form must be 
algebraicly degenerate, with vanishing determinant and therefore matrix rank 2 (since 
an antisymmetric matrix cannot have even rank), a condition that is expressible by the 
algebraic restriction 

ZV^-UJpa] = . (125) 

This has the obvious corollary that the right hand side of (121) must vanish, and hence 
that the helicity current rf will indeed be conserved, i.e. we shall have 

V/ = 0. (126) 

Another, more widely known, consequence of the degeneracy property (125) is that 
(exactly as in the relativistic case explained elsewhere [3]) the vorticity 2-form will 
be orthogonal to a 2 dimensional tangent element, containing the flow vector n M as 
well as the vorticity 3- vector defined by (117), and that (as in the analogous case 
of magnetic 2-surfaces in a perfectly conducting plasma) these 2-dimensional tangent 
elements will be integrable in the sense of meshing together to form well defined vorticity 
flux 2-surfaces. 

The analogue for a 2-form of the formula (85) for Lie derivation with respect 
to a flow field n M takes the form 

nLw^ = ZrfV^w^ - 2V^(w u]a n (T ) , (127) 

in which the first term will drop out identically by the closure property, i.e. the vanishing 
of the exterior derivative 3V[ CT tz7 M „] of the vorticity as an automatic consequence of its 
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exactness property (114). When the field equation (123) is satisfied, it can be seen that 
the second term in (127) will also drop out. We thus obtain the covariant generalisation 
of the well known Kelvin vorticity conservation theorem [19] to the effect that the the 
vorticity 2- form will be conserved by Lie transport, with respect to any arbitrarily 
rescaled multiple of the flow vector, i.e. 

(Cn)£tJ7^ = 0, (128) 

for an arbitrarily variable scalar field (. 

It is to be noted that if, instead of restricting the variation 5n u to be given by 
the worldline displacement formula (86), one merely imposes current conservation by 
adding a Lagrange multiplier term pS/^n^ to the action density, then one will get a 
more restricted dynamical equation to the effect that the momentum covector should 
be the gradient of the Lagrange multiplier ip and thus that it should be irrotational: 

= => uj^ = 0. (129) 

A solution of this irrotational type is the only kind that is allowed in the special case of 
a simple superfluid (on a mesoscopic - i.e. intervortex - scale) for which the scalar p is 
interpretable as being proportional to the quantum phase angle of a bosonic condensate. 

In terms of the Newton-Cartan connection given by (25) the dynamical equation 
(123) can be rewritten in the manifestly gauge invariant form 

u v D v u» = -—Y v V u x ■ (130) 
m 

This is just a covariant reformulation of the well known Euler equation, which is tradi- 
tionally expressed in terms of the pressure function (98), whose variation can be seen 
from (97) to be given in terms of that of the chemical potential x by 

5P = n5x- (131) 
It is thereby possible to rewrite (130) as 

u "V vU i* = - 7 ^(v„0 + -V„P) , (132) 
v P 1 

which can immediately be translated into Aristotelian coordinate notation to give the 
original Eulerian version in the form [19] 

V Ui + tPVjVi = -V,0 - -V t P . (133) 

While this last version has the advantage of familiarity, and (130) has the advantage 
of manifest gauge covariance (with respect to linear Galilean and non-linear Milne 
transformations) it is the version (123) that is most convenient for many mathematical 
purposes, since it involves only exterior differentiation, and can therefore be evaluated 
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in arbitrarily curved (e.g. comoving) coordinates using only partial differentiation, 
without reference to any of the various relevant connections (the oj v p that is covariant 
but curved or the connection that is flat but gauge dependent). An example is 
the demonstration above of the way the use of (123) greatly facilitates the treatment of 
helicity conservation, a concept that is almost trivial (actually simpler than the concept 
of vorticity conservation) in the covariant formalism developed here, but whose original 
derivation, in the frame dependent notation of (133) was of a technical complexity 
such that it was finally obtained (by Moreau and Moffat in the 1960's [15, 16]) about 
a century later than the development of the more elementary precursor concept of 
vorticity (by nineteenth century pionneers such as Kelvin). The advantage of the fully 
covariant approach will be even greater when we go on from simple to multiconstituent 
fluids. 

9. Multiconstituent fluid models 

We now extend the discussion to cases involving several independent - but not 
always independently conserved - currents with current 4-vectors that we shall denote 
by n v where x is a "chemical" index with values ranging over the labels of the various 
constituents involved. In particular the neutron star application for which this work is 
particularly intended, will involve a neutron number density current and a proton 
number density current so in this case the index x will range over the pair of values 
x = n and x = P- Although the total baryon number current n h = n n + n p will be 
conserved, in applications dealing with long term evolution the neutron and proton 
currents will not be separately conserved due to the possibility of transfer of baryons 
from one to the other by weak interactions. To deal with such cases it may be necessary 
to allow for the possibility that a particular current may be characterized by a non 
vanishing value of the destruction rate (per unit spacetime volume) that is defined (as 
the negative of the corresponding creation rate) by 

= -V,< , (134) 

a formula in which, as explained above, it makes no difference what frame may have 
been used to specify the connection involved in the covariant differentiation operator 

The obviously natural way to set up an appropriate Lagrangian for a multicon- 
stituent fluid model is to take a combination of the same general form 

A = A mat +A pot , (135) 

as before, with the material contribution A mat again given by a decomposition of the 
form (105) as a sum of kinetic and internal contributions. As always, the gravitational 
potential energy contribution will simply be given by 

A ot = -0P, (136) 
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where p is the total mass density. However this will now be given as a sum over 
constituents of the form 

p = I]m x n x , (137) 

X 

in which m x is the Newtonian mass per particle associated with the current n£. The 
prescription (136) evidently has the same general form (94) as before when written 
covariantly in terms of the corresponding total mass current vector 

p v =Y.m*n v ^. (138) 

One of the basic principles of Newtonian theory is that although the different contri- 
butions need not be separately conserved (as matter can be transferred from one to 
another by chemical or nuclear reactions) the total mass current (including all relevant 
contributions) will still have to obey the conservation law (108). 

In summation formulae such as (137) and (138) (in which it would be legitimate 
to use the standard shorthand summation convention whereby the explicit use of the 
summation symbol £ is omitted) it is to be noticed that the constituent indices of 
the masses have been written "upstairs" to indicate their contravariant character with 
respect to linear constituent recombinations, in contrast with the currents, with indices 
"downstairs", which undergo recombinations of the corresponding covariant (inverse) 
form. The formulae (137) and (138) can be seen to be covariant, while the resulting 
sums p and p M themselves are actually invariant, when such linear transformations of 
chemical basis are carried out. 

A simple illustration of a change of chemical basis is provided by typical astrophys- 
ical applications for which it may be sufficient to treat the relevant matter (e.g. in a 
stellar atmosphere) as a mixture of hydrogen (with atomic nucleus containing just one 
proton) and helium (with atomic nucleus consisting of 2 protons and 2 neutrons), so that 
in terms of chemical index values x = h and x = He the total mass density will be given 
by p = m H n H +m Hc n Ha . In an equivalent description based on the underlying proton and 
neutron number densities, using index values x = P and x = n, the total mass density 
will be given by an expression of the exactly analogous form p = m p n p +m n n n in which 
the relevant densities are provided by a chemical basis transformation that is given by 
the relations n p = n H + 2n Hc and n n = 2n Hc . Having voluntarily chosen to use down- 
stairs chemical indices for the currents, we have no option but to use upstairs indices 
for the masses because their corresponding transformation will be of contravariant kind 
(i.e. given by the inverse of the covariant transformation matrix) that in this particular 
illustration will be specified by the relations m H = m p , m Ho = 2m p + 2m n . At the cost 
perhaps of obscuring other relevant information, such a change of chemical basis (from 
the atomic reference system to the nuclear reference system) would have the advantage 
of facilitating the exploitation of the empirical fact that for many practical applications 
(particularly in contexts for which a Newtonian description is sufficiently accurate) it 
will be a good enough approximation to take m n ~ m p (a relation attributable to the 
corresponding, but still theoretically unexplained, approximate inequality between up 
and down - but not strange - quark masses). 
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Having dealt with the gravitational potential contribution, we now turn to the ki- 
netic contribution in the decomposition (105) of A kin , for which a multiconstituent 
version can be obtained simply by adding up the contributions, as specified by (96), of 
the separate constituents. We thus obtain a prescription of the form 

A kin = En x X, (139) 

in which, as in (138), the sum on the right is covariant with respect to spacetime 
coordinate transformations. However despite its neat appearance (but due to the non- 
linearity in the defining formula (45) for relevant kinetic momentum covectors p*) this 
contribution (139) is neither gauge invariant (with respect to non-linear Milne or even 
linear Galilean transformations) nor invariant under changes of chemical basis. 

In the manner shown by (70) we can recover gauge invariance of the global inte- 
grated action perturbation (though not of the local unperturbed action density) by 
combining the kinetic contribution with the gravitational potential contribution. This 
gauge invariance at the global perturbation level will evidently be preserved by the 
addition of the extra locally Galilei (and hence a fortiori also Milne) invariant term A. nt 
that is needed to give a chemically invariant value for the total material contribution 
A mat , and hence also for the complete Lagrangian (135). The variation 

M mat = E^<K, ( 14 °) 

of the chemically invariant material Lagrangian will then define a chemically contravari- 
ant set of material momentum covectors yU*. These will be decomposible in the form 

Vl=Pl + Xl, (141) 
where the internal momentum contributions are defined by the variation 

M int = E X *5<, (142) 

of the internal Lagrangian contribution, whose chemical basis dependence endows the 
momentum contributions xt w hh corresponding bad (meaning non contravariant) be- 
havior under chemical basis transformations so as to cancel the bad behavior of the 
kinetic momentum contributions in such a way that the total (141) is chemically 
well behaved. However, although their chemical transformation behaviour is compli- 
cated, the internal momentum contributions have the convenient redeeming feature 
that (unlike the chemically well behaved total /i*) they are automatically invariant with 
respect to Milne (and therefore a fortiori Galilean) gauge transformations, as a conse- 
quence of the postulated gauge invariance of the Lagrangian contribution A int itself. 
In the manner to be derived in a subsequent article [18], this gauge invariance entails 
corresponding Noether identities that are expressible as 

EV^xW^O. (143) 
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and 

En^xl = 0. (144) 

The space projected parts of the internal momentum give rise to the effect (of a 
kind that is familiar in the case of ordinary electron currents in a metallic conductor) 
that is known as "entrainment" whereby the constituent momentum directions may 
deviate from those of the corresponding velocities. However it follows from (143) that 
the deviations will cancel out in the total, so that the (gauge dependent but chemi- 
cally invariant) 3-momentum density, i.e. the space projected part of the 4-momentum 
density (138), as defined by 

W = -f vf r, (145) 
will be expressible as a sum over separate material momentum contributions in the form 



n» = En x », (146) 

X 

in which the individual contributions are given by 

77^ = ^*7^. (147) 

In a manner analogous to that by which the ordinary pressure function P was in- 
troduced by (98) in the simple perfect fluid case, it is useful to define a corresponding 
gauge invariant (and chemically invariant) generalized pressure function ^ for the mul- 
ticonstituent case by the specification 

* = A int -£< X x „, (148) 

It can be seen from (142) that the corresponding infinitesimal variation formula will 
have the form 

5* = (149) 

Putting the kinetic, internal, and external gravitational potential contributions to- 
gether, we now see that the variation of the complete Lagrangian (135) will have the 
(chemically covariant) form 

5 A = EK$K- P 5( t>' ( 15 °) 

while the corresponding variation of the pressure function will be expressible in the 
form 

5V = -EnZ6n*-p6<j>, (151) 

in which the complete (chemically contravariant) particle momentum covectors of the 
various constituent currents are given by expressions of the same form (102) as in the 
single constituent case, namely 

7r- = ^-m x 0V (152) 
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It is to be observed that since the extra (gravitational) term here is purely temporal, 
it makes no difference to the three momentum, so (147) can just as well be written in 
the form 

TJxm = n x 7r^7^ . (153) 

As in the case (106) of a simple perfect fluid, we can rewrite the complete Lagrangian 
in terms of the momenta 7r x and the generalized pressure function \I/ in the form 

A=X>>* + ^. (154) 

Since, as in the single constituent case, the effect of an infinitesimal gauge transforma- 
tion on the momenta will be given in terms of the relevant boost potential (3 by an 
expression of the simple form 

d7r x = -ra x V M /3, (155) 

while the other quantities in (154) will remain invariant, it is evident that the resulting 
infinitesimal gauge variation of the Lagrangian will be given by the formula 

dA = -Em x <Vw3. (156) 

x x 

Subject to the (chemically invariant) restriction that each of the currents should be 
separately conserved, 

V v n£ = 0, (157) 

it can be seen to follow that dA will be expressible as a pure divergence of exactly the 
same simple form (109) as in the single constituent case. Since it can be seen from (87) 
that the restriction (157) is preserved by the generic flow displacement (84) it follows 
that the corresponding variational equations of motion that we shall derive below will 
be gauge independent. 

It is to be observed that the condition (109) for the gauge invariance of the integrated 
action will still be satisfied even if the currents do not satisfy the separate conservation 
conditions (157) but are restricted only by the single conservation condition (108), that 
must always be satisfied by the total Newtonian mass current (138). However, if the 
currents are not separately conserved, the total mass conservation condition (108) will 
not automatically be preserved by a generic set of independent flow displacements of 
the form (84), so the variational principle will no longer provide an automatically gauge 
invariant set of field equations. 

When each current is subject to its own independent displacement the gen- 
eralization of (110) that we finally obtain by substituting (86) in (150) will take the 
form 

SA = V M (2Evr>[^ ] ) - - PH, (158) 

in which, for each constituent, the covector / x is interpretable as the 4- force density 
acting on the corresponding current n£, and in which the value of this force density 
can be read out as 

/ x = 271^7^] + 7T X M V,< . (159) 
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It is this that must vanish in the strictly conservative case for which the variational 
field equations are satisfied. 

Whenever the separate conservation conditions (157) actually are satisfied, the force 
density will reduce to the simple form 

f£ = <w\,, (160) 

using the notation 

= SV^] , (161) 

for the generalized vorticity tensors, which unlike the momenta from which they are 
derived, can be seen from the form of (155) to be gauge invariant, 

dwV = 0. (162) 

Although the decay rates (134) are also gauge invariant, 

d£> x = 0, (163) 

whenever they are non zero, i.e. when the separate conservation conditions (157) are 
not satisfied, the resulting force density contributions 

/* = n>V - n^V x , (164) 

will no longer be gauge invariant, but will transform according to the rule 

df x =m x P x V^. (165) 

It can be seen from the general formula (127) that when the separate current conser- 
vation conditions (157) are satisfied, the Lie derivatives of the vorticities with respect to 
the corresponding flow fields will be given in terms of the corresponding force densities 
by 

n x £^V = 2V [M / X ] . (166) 

When the full set of variational field equations is satisfied, so that the forces densities 
/* all vanish, it can be seen that each vorticity will be conserved by transport along the 
corresponding flow lines, i.e. each constituent will satisfy a Kelvin type conservation 
law of the form (128), namely 

(C x n x )£^V = 0. (167) 

for arbitrary scalar fields £ x - 

As in the corresponding relativistic case [10], one can go on to generalize the single 
constituent helicity vector (119) to a vector valued helicity matrix defined by 

f Y " = i^vr^V . (168) 
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The antisymmetric part of this matrix will be exact in the sense of having the form of 
a divergence, 

^ = ' Vl/ (^^) , (169) 

of an antisymmetric tensor, with the implication that it will automatically be closed in 
the sense that its own divergence will vanish identically, i.e. 

v ^ [xy] m = o . (170) 

It follows that if a constituent with label x say is characterized by the property of 
irrotationality, meaning that tu x ^ u vanishes, and more particularly, as in (129), in the 
case of a constituent that is superfluid, and thus characterized (on a mesoscopic scale) 
by a momentum covector that is the gradient of a corresponding phase scalar </? x , so 
that for any other label value y say the corresponding helicity matrix component r] YXtl 
will also vanish, i.e. ^ YXAt = 0, then the corresponding transposed component will 
automatically be conserved, i.e. we shall have 

^ = V^ X wV = ° =>• V M 7f YM = 0. (171) 

Regardless of any irrotationality constraint that may be satisfied, it can be seen - 
for the same reason as in the single constituent case characterized by (126) - that the 
vanishing of the force density covector in (160) will always be sufficient to ensure 
that the divergence of the corresponding diagonal component of the helicity matrix 
will vanish, i.e. that the constituent under consideration will be subject to a helicity 
conservation law [2, 10] having the form 

V M rf^ = 0. (172) 
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Appendix: helicity current in the traditional 3+1 formalism 

The efficacity of 4-dimensionally covariant treatment is clearly demonstrated by the 
ease with which the foregoing helicity conservation laws have been obtained as almost 
obvious consequences of the dynamical equation. This contrasts with treatment using a 
3+1 space versus time decomposition, in which an equivalent derivation of the helicity 
conservation law [16] requires much greater algebraic effort and ingenuity, as shown in 
this appendix. 

The time and space components of the helicity current are found from (119) and 



the specification (6) to be expressible as 

V° = -e^TTiVjiTk , (173) 

jf = Tr e ijk VjTr k - s^tt^tt, + e^%V k n Q . (174) 

In the following - as is usual within an Aristotelian - Cartesian framework. - the 



space indices will be replaced by the familiar arrow notation. We also introduce the 
cross-product ((f) x tp) 1 between two forms fa and <pi as the contravariant vector e tjk (f)j(p k . 
The curl and 3-divergence of a vector V are defined respectively as (curl V) % = e tjk VjV k 
and div V = ViV\ The time derivative will be written as V = d t . 
With this notation, the helicity 4-vector will be given by 

T]° = — tt .curl 7T , (175) 

ff = 7T CUrl TT — TT X d t TT + TT X VtT () . (176) 

Considering the example of the single perfect fluid model, the components of the 
4-momentum covector (102) are tt q = —S and tt = mv, in which we have introduced 
the total particle energy £ = \mv 2 + mcf) + \- Hence the helicity current reduces to 

rf = -2m 2 Q.v, (177) 

ff = —2m£uj + V£ x mv + m 2 d t v x v (178) 

— # 

in terms of the kinematic local angular velocity u = Neurit? (divu; = 0). 

The conservation law of this helicity current will now be derived within the 3+1 
spacetime decomposition, starting from Euler's equation in the form [19] 

md t v + V£ + 2mujxv = 0, (179) 

whose contraction with mus yields 

m 2 u.d t v + mu.V£ = 0. (180) 

Taking the curl of the Euler equation followed by the dot product with the velocity 
gives 

m 2 v.d t uj + m 2 v.cni\{Q x v) = . (181) 
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Adding the two equalities and simplifying, we obtain 

2m 2 d t (u.v) + 2div(mu£) + 2m 2 (v.u)div v + 2m 2 v.(v.Vu) - 2m 2 v.(Q.Wv) = . (182) 

Taking the cross-product of the Euler equation with the velocity, the divergence and 
using the standard identity 

ViW.V) = W x curl V + V x curl W + V.VW + W.VV , (183) 

we obtain 

div(m 2 v x d t v + mv x V£) + 8m 2 u.(t? x u) + Am 2 Co .{v .V v) 

-2m 2 (v.Lu)divv-2m 2 v.(v x curld;) - 2m 2 v.(vWCu) - 2m 2 v.(u.Vv) = . (184) 

Combining this with (182), one obtains the space-time decomposition of (126), namely 
the local version of the original Moreau-Moffat helicity conservation law [15, 16], in the 
form 

2m 2 d t (u).v) = dlv(-2muj£ + m 2 d t v x v + m(V£) x v^j . (185) 
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